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Abstract 

We report new aspects of the recent theoretical progress in heavy quarkonium physics. 
(1) Contrary to wide beliefs, the gross structure of the bottomonium spectrum is described 
well by the non-relativistic boundstate theory based on perturbative QCD. (2) This leads 
to a new physical picture of the bottomonium states: the boundstate mass is composed 
mainly of the self-energies of b and b accumulated inside the boundstate. (3) A connection 
to the conventional phenomenological potential-model approaches is provided. 



*Talk given at "Accelerator and Particle Physics Institute (APPI 2001)", Morioka, Japan, Feb. 20-22, 2001. 



1 Introduction 



In this paper we review the results of the recent studies on the spectra of heavy quarkonia 
(bottomonium, charmonium and B c states) Jl|, |2|. We focus mainly on the bottomonium 
spectrum below; the results of the charmonium and B c spectra will be summarized at the end 
of the paper. 

For over 20 years, major theoretical approaches to the charmonium and bottomonium spec- 
troscopy have been those based on various phenomenological potential models. In each of the 
models a non-relativistic Hamiltonian, 

# = |^ + Veno(r), (1) 

is assumed and a phenomenological potential lpheno( r ) is determined such that the observed 
quarkonium spectra (and some other physical observables) are reproduced, starting from a 
simple ansatz for the form of the potential. The determined potentials in all models have 
more or less similar slopes in the range 0.5 GeV" 1 < r < 5 GeV" 1 , which may be represented 
by a logarithmic potential oc logr + const. These phenomenological-model approaches have 
successfully elucidated nature of the quarkonium systems, such as their leptonic widths and 
transitions among different levels, besides reproducing the energy levels. See e.g. Ref. |J for a 
most recent analysis based on the potential models. An apparent deficit of these approaches is, 
however, a difficulty in relating phenomenological parameters to the fundamental parameters 
of QCD. 

The reason why people had to resort to phenomenological models is because the theory of 
non-relativistic boundstates, which has been successful in describing the spectra of the QED 
boundstates, failed to reproduce the charmonium and bottomonium spectra in QCD. Within 
this theory, the quarkonium states in the leading approximation are described by the Hamilto- 
nian of non-relativistic quantum mechanics with a Coulomb potential 

V c (r) = -C F ^. (2) 
r 

It stems from one gluon exchange between the quark and antiquark, where CV = 4/3isa color 
factor. Compare the observed bottomonium spectrum and the Coulomb spectrum shown in 
Fig. |l](a). In the Coulomb spectrum, the level spacing between consecutive nS states decreases 
rapidly as 1/n 2 . On the other hand in the bottomonium spectrum the level spacing appears 
roughly constant; furthermore, the separations between the S and P states as well as the fine 
splittings among nPj states are sizable. Thus, the level structures look qualitatively very differ- 
ent. Including higher order corrections, the Coulomb potential changes to the QCD potential, 
which is given roughly by replacing as in the numerator of Eq. (D) with the running coupling 
constant evaluated at scale 1/r, i.e. as (yU = 1/V). Accordingly the QCD potential is bent down- 
wards at long distances as compared to the Coulomb potential. As can be seen in Fig. |l|(b), 
the QCD potential becomes singular at a fairly short distance. At r > 0.2 GeV -1 , the QCD 
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(a) (b) 

Fi gure 1: (a) A comparison of the Coulomb spectrum and the observed bottomonium spectrum. The pa- 
rameters nib and as of the Coulomb spectrum are adjusted to reproduce the observed 15* and 2S levels, (b) 
A comparison of the renormalization-group-improved QCD potential, the Coulomb potential and a typical 
phenomcnological potential. 



potential is poorly convergent as the higher order corrections are included. There seems to be 
no chance that it explains the phenomenologically determined potentials, which deviate from 
the Coulomb potential in the upward (i.e. opposite) direction. The large discrepancy of the 
QCD potential and phenomenological potentials, as well as the poor convergence of the QCD 
potential, have been thought as indications of large non-perturbative effects inherent in the 
heavy quarkonium systems. In fact the difference between a typical phenomenological poten- 
tial and the Coulomb potential tends to be a linearly rising potential at distances r > 1 GeV -1 , 
suggesting confinement of quarks. 

During the last few years, the theory of non-relativistic QCD boundstates has developed re- 
markably. There were two important developments: (1) The complete next-to-next-to-leading 
order corrections to the energy levels have been computed [§], §, || [7|, f|. Also, some of the 
(non-trivial) next-to-next-to-next-to- leading order corrections have been calculated || [lO], |TT|1 . 



(2) The renormalon cancellation in the energy levels was discovered fT2| , |13fl . The upshot is 
that convergence of the perturbative expansions of the energy levels improves drastically if we 
express the levels in terms of the MS mass instead of the pole mass of a quark. These theoreti- 
cal developments enabled accurate perturbative computations of the quarkonium energy levels. 
Let us demonstrate the improvement of convergence for the T(IS') and T(2S) states: 
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T(15): For /i = 2.49 GeV, a 3 (jJ,) = 0.274, mf s (mf s ) = 4.20 GeV/m bjPolc = 4.97 GeV, 

M T (is) = 9.94- 0.17- 0.20 - 0.30 GeV (Pole-mass scheme) (3) 
= 8.41 + 0.84 + 0.20 + 0.013 GeV (MS-scheme) . (4) 



• T{2S): For li = 1.09 GeV, a s (/x) = 0.433, mf s {mf s ) = 4.20 GeV/m bjPole = 4.97 GeV, 

M T(25) = 9.94- 0.10 - 0.19 -0.45 GeV (Pole-mass scheme) (5) 
= 8.41 + 1.46 + 0.093 + 0.009 GeV (MS-scheme). (6) 

As can be seen, when the pole mass is used the series are not converging, whereas the series 
show healthy convergent behaviors when the MS mass is used. As important applications of 
these developments up to date, the theory enabled precise determinations of the MS-mass of 
the bottom quark 0, [14], [| [L5| and (in the future) of the top quark (16| from (mainly) the 



energy levels of the lowest-lying states. The main uncertainty comes, in the bottomonium case, 
from the (essentially) unknown non-perturbative contributions. These are generally claimed to 
be around 100 MeV and ultimately set the precision of the prediction. 

Based on the above theoretical developments we analyze the consistency between the whole 
level structure of the quarkonium as predicted by the boundstate theory (perturbative QCD) 
and that of the experimental data. Now that we can make accurate predictions, we can ex- 
tract upper bounds on the non-perturbative contributions to the energy levels by comparing 
the perturbative predictions, at the current best accuracy, with the experimental data. It 
turns out that non-perturbative effects should be much smaller than what have been believed 
conventionally. 

The paper is organized as follows. In Sec. 2 we review briefly the necessary theoretical 
framework. We analyze the bottomonium spectrum numerically in Sec. 3 and estimate the 
errors of the theoretical predictions in Sec. 4. We interpret the result and provide a new 
physical picture of the bottomonium states in Sec. 5. Then we discuss a connection between our 
approach and the phenomeno logical potential- model approaches in Sec. 6. We draw conclusions 
in Sec. 7. 



2 Theoretical Framework 

2.1 1/c expansion 

We state briefly the theoretical framework used in contemporary calculations of the spectrum 
of heavy quarkonia. The problem is reduced to a quantum mechanical one as follows. We first 
compute the quark- ant iquark (off-shell) scattering amplitude in ordinary perturbative QCD. We 
then determine a quantum mechanical Hamiltonian such that the quark-antiquark scattering 
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amplitude computed within quantum mechanics matches the former amplitude order by order 
in expansion in 1/c (inverse of the speed of light): 



H = H + - Hi + \ H 2 + ■ ■ ■ ■ (7) 



v c 2 



The expansion is a double expansion in as = g s /(4:Tchc) and (3 = v/c. Then we solve the 
non-relativistic Schrodinger equation 

Hip n (r) = E n ip n (r) (8) 

to determine the boundstate wave functions and energy spectrum, order by order in 1/c expan- 
sion. Provided the quark and antiquark inside the heavy quarkonium are non-relativistic, the 
expansion in 1/c leads to a reasonable systematic approximation. 

It should be noted that this procedure parallels a more familiar procedure for the calculation 
of the masses of one-particle states. The matrix element of the above Hamiltonian (p\H \p') 
is, in the language of perturbative QCD, the sum of two-particle irreducible diagrams.^ One 
may then compare Eq. (H) with e.g. the computation of the photon and Z boson masses by 
solving the eigenvalue equation 

_ / n 77 ( s ) n 7Z ( s ) \ _ 

^ Yi Zl {s) Yi zz {s) ) ~ u ' 

where Ily(s) represent the sums of one-particle irreducible diagrams [U zz (s) includes M ( z ee) }. 
Presently the Hamiltonian is known up to 0(l/c 2 ) in the Coulomb gauge || 18] : 



H = ^-C F ^, (9) 
m r 

Hi = -C F ^-(g)-{/3 log( / i'V) + a 1 }, (10) 

H 2 = -J^-^V'© 2 ^•{/ 5 o[log 2 (/i , V) + ^] + (/3 1 + 2/? a 1 )log(/i'V) + a 2 } 
7rC F a s 3 3C F a s =» - C F a s f ^ 2 , 1 \ C A C F a\ 



m? 2m 2 r 3 2m 2 r \ r 2 % 3 ' 3 1 ) 2mr 5 

C F a s fS 2 JS-f) 2 47r, oc , 2 & 



2mMr3 V-T^-W}. (») 

where m denotes the pole mass of the quark; as = (Xsip)] C F = 4/3, Ca = 3 are color factors; 
fi' = // e lE . The lowest-order Hamiltonian Ho is nothing but that of two equal- mass particles 
interacting via the Coulomb potential. 



This is true with respect to the diagrams in the time-ordered (old-fashioned) perturbation theory pj 
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(a) 



(b) 



Figure 2: The QCD potential in the large-/?o approximation truncated at 0(ctg +1 ) term. We set \x = 2.49 GeV, 

rii = 4 and ots{n) = 0.273 [corresponding to af (Mz) = 0.1181]. (a) Before subtraction of the leading 
renormalon. (b) After subtraction of the leading renormalon. These figures are taken from 

2.2 Renormalon cancellation 

In addition to the first three terms Hq, Hi, H2 of the Hamiltonian, part of the higher order 
terms H n are known. Specifically, these are part of the static QCD potential Vqcd(?"). From 
their analysis, it has been known [[HJ that the series expansion of Vqcd(^) in ®-s{n) diverges 
rapidly at high orders, and that this results in an uncertainty of Vqcb(t) of order Aq CD even 
within perturbative QCD. This problem is referred to as the renormalon problem. This is 
demonstrated in Fig. 0(a), where the QCD potential in the so-called "large-/3o approximation" 
is shown up to the (JV+l)-th term, J2n=o V8o\ r )i f° r N = 0,1,2, .. . and n\ = 4. We see that 
the higher order corrections are indeed large and almost constant (independent of r). 

It was found fi"2| , |i~3|1 that the leading renormalons contained in the pole mass |2(J and in 
the QCD potential |L9| cancel in the total energy of a static quark- ant iquark pair if the pole 
mass m po ie is expressed in terms of the MS mass: 



where q = \q\. The above equations show that the potential Vqcd(^) is essentially the Fourier 
transform of the Coulomb gluon propagator exchanged between quark and antiquark, and that 
the difference of m po i e and is essentially the infrared portion of the quark self-energy. The 
renormalon contributions originate from the infrared region, q ~ Aqcd ; of the loop integrations, 



E tot (r) = 2m pole + VqcdM 



(12) 




(13) 



(14) 
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where the running coupling constant ats{q) becomes large. Namely the infrared gluons cause bad 
behaviors of the perturbation series at high orders. The signs of the renormalon contributions 
are opposite between VqcdO") anc ^ m poie because the color charges are opposite between quark 
and antiquark while the self-enregy is proportional to the square of a same charge. Their 
magnitudes differ by a factor of two because both the quark and antiquark propagator poles 
contribute in the calculation of the potential whereas only one of the two contributes in the 
calculation of the self-energy. Expanding the Fourier factor e iq ' r in a Taylor series for small q, 
the leading renormalon contributions cancel in E tot (r). 

As a result of this cancellation, the series expansion of the total energy in as(fj) converges 
better if we use the MS mass instead of the pole mass. In order to demonstrate the improvement 
of convergence, we show in Fig. |2](b) the QCD potential in the large-/5o approximation [up to 
the (iV + l)-th term] after the leading renormalon is subtracted at each order of «s(/i). One 
sees that the series expansion of the potential has become much more convergent as compared 
to Fig. ^(a). Note that the higher-order corrections raise the potential at long distances after 
subtraction of the renormalon contribution. Some important aspects are: 

• The pole mass of a quark is ill-defined beyond perturbation theory. It can be determined 
only when the quark can propagate an infinite distance. Generally accepted belief is 
that when quark and antiquark are separated beyond a distance ~ Aqq D the color flux 
is spanned between the two charges due to non-perturbative effects and the free quark 
picture is no longer valid. On the other hand, the total energy (or the mass) of a quarko- 
nium, which is a color-singlet state, is physically meaningful. A color-singlet state can 
propagate for a long time and the notion of its mass is not limited by the hadronization 
scale. 

• When the size of a color-singlet system is much smaller than Aq^ d , infrared gluons with 
wavelengths Aq^ d cannot couple to color sources inside the system — such a picture is 
naturally described by the bare QCD Lagrangian. Hence, if we use ttlj^, which is more 
closely related to the bare mass than m po i e is, contributions from the infrared gluons 
vanish in E tot (r). 



See e.g. ET] for an introductory review of the renormalons in the heavy quarkonium states. 



3 Bottomonium Spectrum: A Numerical Analysis 

In this section we examine the bottomonium spectra numerically Jl| . According to the formal- 
ism explained in the previous section, the energy levels are computed analytically as functions 
of Qg(/i), fi and mf, = m^ /[S (m^ IS ) (the 6-quark MS mass renormalized at the MS-mass scale). 
The dependence on the scale /i arises from truncation of the series expansion at finite order. 
In this section we set the number of massless flavors as rii = 4 and neglect the effects of the 
non-zero charm quark mass in the bb systems. 

The algorithm of our calculations goes as follows. 
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1. We take the strong coupling constant with the present world average value as input ,[] 
af{M z ) = 0.1181 ± 0.0020 J22 . 



2. The scale \i is fixed for each state X from the minimal sensitivity condition: 

= 0. (15) 

3. We fix nib from the mass of the vector ground state: 

Er(is)(t*is, as(n is ),m b ) = E^ 1S) = 9.460 GeV. (16) 

A comment regarding the scale determined by the scale-fixing prescription Eq. ( |T5| ) is in 
order. We find that, for most of the bottomonium states, the convergence properties of the 
series become optimal at fi ~ fix, and that the scale becomes close to the inverse of the physical 
size of the boundstate X. If the scale fixed by Eq. (|15D evidently does not fulfill these conditions, 
then the theoretical predictions obtained in this way will be considered unreliable. We will show 
that this typically happens for the higher levels, where the coupling constant becomes bigger 
than one. 

The 6-quark MS mass fixed by Eq. ([16]) is given by 



af{M z ) 


0.1161 


0.1181 


0.1201 


m b 


4.221 GeV 


4.203 GeV 


4.184 GeV 



Using these masses as inputs and Eq. (|15|), we can calculate the energy levels of other observed 
quarkonium states. These are shown in Figs. []. Only those levels which can be predicted reliably 
are displayed. The level spacings become wider for larger otg\Mz), which is consistent with our 
naive expectation. If we take an average of the 5-wave and P-wave levels corresponding to each 
principal quantum number n, the theoretical predictions with a^(Mz) = 0.1181 reproduce the 
experimental values fairly well. On the other hand, the predictions for the S—P splittings and 
the fine splittings are smaller than the experimental values. 

Also the numerical values of the series expansions of the energy levels are shown in Tab. [l] 
in the case o§\m z ) = 0.1181. In this case, the series expansions for the l 3 Si, l 3 Pj, 2 3 S' 1 , 2 3 P 
and 3 3 S"i bottomonium states converge well. For these states the differences of the theoretical 
predictions and the experimental data are typically 30-60 MeV. Convergence of the series 
expansions is poor for the 2Pi, 2P2 and 4S states. We consider that the theoretical predictions 
for these levels are unreliable and marked the corresponding theoretical predictions with sharps 
(Jj). The differences E e x P — Ex are rather large for the states with sharps. Notice that for 
these states the corresponding as(^x) becomes larger than one, indicating a breakdown of 

' We evolve the coupling and match it to the coupling of the theory with n; = 4 via 4-loop running. 



^-E x {^as{^),m h ) 
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Coulomb 



Experimental 
Values 



Perturbative QCD 
a s (M z ) = 0.1181 



Experimental Perturbative QCD Experimental Perturbative QCD 

ValueS a s (M z ) = 0.1161 Values a s (M z ) = 0.1201 



Fi gure 3: Comparisons of the theoretical predictions and the experimental data for the bottomonium spectrum. 
The input values of the theoretical predictions are a^(Mz) = 0.1181, 0.1161 and 0.1201. We set m c = and 
the number of massless flavors as n/ = 4. The solid and dashed lines represent the S- and P-states, respectively. 
We show only those levels which we can compute reliably. 
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State X 


rpexp 

&x 


Ex 


E x p -E x 






7+3 
a x 






T(l 3 S'i) 


9.460 


9.460 





0.837 


0.204 


0.013 


2.49 


0.274 


T(PP ) 


r\ Odd 

9.860 


9.905 


—0.045 


1.38 


0.115 


0.003 


1.18 


C\ A C\C\ 

0.409 


T(l 3 Pl) 

T(1 3 P 2 ) 


9.893 


9.904 


-0.011 


1.40 


0.098 


0.002 


1.15 


0.416 


9.913 


9.916 


-0.003 


1.42 


0.086 


0.003 


1.13 


0.422 


T(2 3 Si) 
T(2 3 P ) 
T(2 3 Pi) 
T(2 3 P 2 ) 
T(3 3 ^) 


10.023 


9.966 


+0.057 


1.46 


0.093 


0.009 


1.09 


0.433 


10.232 


10.268 


-0.036 


2.37 


-0.66 


0.15 


0.693 


0.691 


10.255 


10.316 s 


-0.061' 


3.97 


-3.56 


1.50 


0.552 s 


1.20 


10.268 


10.45 7 s 


-0.189 s 


4.55 


-5.03 


2.53 


0.537 s 


1.39 


10.355 


10.327 


+0.028 


2.34 


-0.583 


0.163 


0.698 


0.684 


T(A 3 S 1 ) 


10.580 


11.760 s 


-1.180 s 


5.45 


-6.47 


4.38 


0.527 s 


1.61 



Table 1: Comparisons of the theoretical predictions of perturbative QCD and the experimental data for 
as(Mz) = 0.1181. All dimensionful numbers are in GeV unit. 



the perturbative series. Generally, for states, which we consider reliably calculable, the scale 
dependence decreases as we include more terms of the perturbative series. For states, whose 
predictions we consider unreliable, the series becomes much more convergent if we would choose 
a scale different from (typically larger than) fix- 



4 Error Estimates 

Besides non-perturbative corrections, there are two kinds of uncertainties in our theoretical 
predictions for the bottomonium spectra. These are listed below and in Tab. |2|. Also these 
examinations indicate that the theoretical predictions for the 2 3 Pi, 2 3 P 2 and 4 3 5*i bottomonium 
states are quite unstable, while the prediction for the 2 3 Po state is at the boundary. 

1) Charm Mass Effects. We have also made an analysis of the bottomonium spectrum 
including finite charm mass effects. We will report the full results in a separate paper 



[23fl. Here we only summarize some of the qualitative features of the effects and include 
them as a part of the uncertainties of our present analysis. The corrections to the reliable 
predictions turn out to be positive and to become larger for higher states, ranging up 
to ~ 200 MeV. Much of the effects, however, can be cancelled by decreasing the input 



a 



(5) 



(Mz) within its present uncertainty.^ 



Uncertainties from Higher- Order Corrections. We take the maximum value of the follow- 
ing five estimates as an estimate of uncertainties from unknown higher-order corrections 
for each series expansion: (i) The difference between the theoretical predictions computed 



* We note that the sensitivities of the higher levels to a variation of a^(Mz) increase by the charm mass 
effects due to the decoupling of the charm quark. 
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Estimates of higher-order corrections 





(i) 


(ii) 


(iii) 


(iv) 


(v) 


±max 


fj, = 2^ x 




—2 


+1 


U 


1 *7 

±7 


+5 


1 *7 

±7 


+ 16 




T(1 3 P ) 


+ 15 


-14 





±3 




±15 


-53 




T(l 3 Pi) 
T(1 3 P 2 ) 
T(2 3 Si) 
T(2 3 P ) 
T(2 3 Pi) 
T(2 3 P 2 ) 
T(3 3 5i) 


+22 


-8 


+7 


±2 




±22 


-48 




+16 


-16 





±3 




±16 


-54 




+ 18 


-19 


-1 


±9 




±19 


-73 




+33 
-4 


-57 
+1637* 


-21 
-66« 


±150 
±1500* 




±150 
±1637" 


-120 
-97 




-136 


+2136 tt 


-206« 


±2530* 




±2530" 


-229 




+37 


-63 


-36 


±163 




±163 


-161 




T(4 3 ,S 1 ) 


-639» 


+2936» 


-1425 tt 


±4380« 




±4380* 


-1361 



Table 2: Variations of the theoretical predictions of Tab. [l] when the uncertainties 2) discussed in Sec. || are 
separately taken into account. All dimensionful numbers are in MeV unit. Those values corresponding to the 
unreliable theoretical predictions are marked with sharps. The input parameter is taken as as(Mz) — 0.1181. 
The column "±max" lists ±max{|(i)|, |(ii)|, |(hi)|, |(iv)|, |(v)|}. In the last column we write conservative 
estimates with the scale choice fj, = 2 fix- 



using as{fi) as obtained by solving the renormalization-group equation perturbatively at 
4 loops (Eqs. (3) and (11) of Ref. [p4[ ) and numerically at 4 loops (the data of Tab. [I]). 
(Note that the number of energy levels that can be determined in a reliable way is larger 
with the former definition of the running coupling constant. Also in that case, reliable pre- 
dictions turn out to be close to the experimental values.) (ii) The difference between the 
theoretical predictions computed using the 3-loop and the 4- loop (as in Tab. [[]) running 
coupling constants, fixing ctg (Mz) = 0.1181. (iii) The difference between the theoreti- 
cal estimates obtained by fixing fix on the minimum of \E^\ and the results of Tab. 0, 
obtained by fixing nx via the condition ([151) . (iv) The contribution ±|£'^' ) | from Tab. |l[ 
(v) For the IS states we consider the 0(a s m) corrections calculated in the large-/9 ap- 



proximation in [11 



For comparison, we list more conservative error estimates. These are the variations of mj 
and Ex when we fix the scale as twice| of the minimal sensitivity scale: \x = 2/ix, where 
fix is determined from Eq. (|T5|). 



§ If we fix the scale as half of the minimal sensitivity scale, fi = /Ux/2, the predictions for the n = 2 
bottomonium states appear to be meaningless, since the scales are quite close to the infrared singularity of the 
runnning coupling constant, and the predictions for the n > 3 states do not exist, since the scales lie below the 
infrared singularity. 
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5 Interpretation 



The most non-trivial feature of the present theoretical predictions for the bottomonium spec- 
trum is that the level spacings between consecutive n's are almost constant, whereas in the 
Coulomb spectrum the level spacings decrease as 1/n 2 . Conventionally, this same difference 
between the Coulomb spectrum and the observed quarkonium spectra motivated people to 
construct various potential models. It is, therefore, imperative to elucidate how the above 
perturbative QCD calculation is able to reproduce such a level structure. We will focus on two 
points: (1) the leading renormalon cancellation, which implies that infrared physics decouples; 
this is essential to obtain convergent series expansions; (2) the meaning of the scale fix chosen 
by the scale-fixing prescription (|T5|). 

According to the discussion in Sec. 2.2, physical meaning of the renormalon cancellation 
can be understood as the decoupling of infrared gluons in the computation of the quarkonium 
energy levels: The gluons, whose wavelengths are much larger than the size of the color-singlet 
boundstate, cannot couple to it. In Eqs. ( |13D and (|14D , the integrands essentially cancel each 
other in the region q < 1/r. We may replace r by the size ax of a boundstate X in Eqs. (|I2|)- 
(14) and write 

^ f d 3 q _ Anotsiq) ,^„ s 

^»+/ p^ c '-?^- < 17 > 

l/ax ^ q<fn 

This expression may be made somewhat more accurate. Let us approximate the energy level 
as 

E x = 2m pole + E hi n,x, (18) 

f d 3 q „ 4:iras(q) 2C F f w , , . 

2m pole ~ 2m + / C F = 2m + — - / dqa s (q), 19 

J {ZTCy q Z 7T JO 

q<m 

E hin ,x ^ (X | + Vqcd I X ) . (20) 

^polc 

Here, Vqcd(<?) — —C F ATias(q)/q 2 is the QCD static potential in momentum space; | X ) denotes 



2Ct? r°° n 2 
Ex ~ 2m+—^ dqas ( q )f x ( q ) + (x\J—\X) (21) 

7T JO m D0 le 



the Coulomb wave function (the zeroth-order approximation). From Eqs. (|l~9|) and fl2~0|) we 
obtain 

2C F foo ^ ^ f f _ s | /vl JP 

^pole 

- 2m+^/ dqa s (q)fx{q). (22) 

71 JO 

The last approximation follows from the fact that the kinetic energy contribution to the bot- 
tomonium levels turns out to be numerically smallQ (notice that this does not contradict the 

*It is about 17% of E { x ] for the X = IS state, 6% of for the X = 2S state, 4% of E^ for the X = 35 
state. Moreover, these contributions tend to cancel each other in the level spacings. 
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1 . 5 



1 



0.5 



0.5 1 1.5 2 2.5 3 
q [GeV] 

Figure 4: The support functions fx(l) vs. q for X = IS 1 , 2S and 3S* (solid lines). fx{l) is calculated using 
m poic — 5 GeV and a different a>s(fJ>x)i taken from Tab. [l], for each X . Vertical lines represent the corresponding 
scales /ix taken from the same table. Also a^^q) is shown by a dashed line. 

virial theorem, since the static potential we are considering here incorporates the running of 
as and therefore is not simply of the form 1/r). The support function fx is given by 

fx{q) = 0{m-q)- f°° drr 2 \R x (r)\ 2 (23) 
Jo qr 

where Rx{f) is the radial part of the Coulomb wave-function of X. fx{q) is almost unity in 
the region l/a x ~ q < 7n, where a x is the inverse of the dumping scale of fx and may be 
interpreted as the size of the boundstate X. Thus, Eq. (|22|) can be identified with Eq. ([[7]) 
qualitatively. 

For the IS state fx(q) dumps slowly as q decreases. For other states fx(q) dumps rapidly 
from scales which are somewhat smaller than the naive expectations (C Ftts^poie) / W • 111 Fig- H 
we show fx for different states calculated with m po i e = 5 GeV and asil^x) taken from Tab. [L| 
The corresponding values of fix are also displayed. For those states which we consider the 
predictions reliable, \ix is located within the range where fx(q) — 1 (close to the infrared 
edge); for those states with unreliable predictions, nx is located out of this range but far in the 
infrared region. This comparison shows a relation between the scale fix and the wave function 
(consequently the size) of the corresponding boundstate X. 

Eq. (P) or (H) tells that the major contribution to the energy levels comes from the region 
l/ax ~ q ~ fn of the self-energy corrections of quark and antiquark (apart from the constant 
contribution 2m). That is, with respect to only the contributions of the gluons inside the 
boundstate, the self-energies of quark and antiquark dominate over the potential energy between 
the two particles. Hence, we find a qualitative picture on the composition of the energy of a 
bottomonium state: 

(I) The energy levels of bottomonium are mainly determined by (i) the MS masses of b and 
b, and (ii) contributions to the self-energies of b and b from gluons with wavelengths 
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Figure 5: The total energy of a heavy quarkonium state is carried by the MS masses of quark and antiquark 
and by the gluons inside the boundstate. In the latter contributions the self-energies of quark and antiquark 
dominate over the potential energy between the two particles. 

1/rn < A < ax- The latter contributions may be regarded as the difference between the 
(state-dependent) constituent quark masses and the current quark masses. 

This picture is reminiscent of the (long-established) interpretation that the masses of light 
hadrons consist of the constituent quark masses. See Fig. |5]. 

Now we are in a position to see why the level spacings among the bottomonium excited 
states are much wider than those of the Coulomb spectrum. We know that if the coupling 
as(q) were independent of q, the energy spectrum would become the Coulomb spectrum. In 
Fig. ^ also ces(q) is shown. We see that as nx increases from 1 to 3, the coupling ces(q), close 
to the dumping scale of fx{q), grows rapidly. According to Eq. fll7| ) or (p2|), as the integral 
region extends down to smaller q, the self-energy contributions grow rapidly in comparison to 
the non-running case. They push up the energy levels of the excited states considerably and 
widen the level spacings among the excited states. Thus, we may draw the following qualitative 
picture of the bottomonium level structure: 

(II) Level spacing between consecutive n's increases rapidly with n as compared to the Coulomb 
spectra. This is because the self-energy contributions (from 1/fn < A < ax) grow rapidly 
as the physical size ax of the boundstate increases. 

6 A Link to Phenomenological Potentials 

We return to the problem which we discussed in the introduction (Sec. 1): How can we un- 
derstand the "linear potential" in the difference between the Coulomb potential and a typical 
phenomenological potential? 

According to the interpretation given in the previous section, the total energy Etot (?") = 
2m po i e + Vqcd(^) of a bb system determines the bulk of the bottomonium spectrum. Hence, we 
examine the series expansion in a,g(/i) of the total energy, expressed in terms of the 6-quark MS 
mass, E tot (r;fnb,as(f J ')), and compare it with phenomenological potentials 0. The obtained 
total energy depends on the scale /i due to truncation of the series at C(a|). One finds that, 
when r is small, the series converges better and the value of E tot (r) is less /^-dependent if we 
choose a large scale for /i, whereas when r is larger, the series converges better and the value of 
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Figure 6: A comparison of the total energy of a bb system and typical phenomenological potentials (solid 
lines). Effects of the non-zero charm quark mass are included. For a reference, we show typical sizes of the 
bottomonium and charmonium S states as determined from the r.m.s. interquark distances with respect to the 
Cornell potential: J (r 2 )„ ... 

1 V V / Cornell 



Etot{f ) is less /i-dependent if we choose a smaller scale for fi. Taking into account this property, 
we will fix the scale \x in two different ways below: 



1. We fix the scale fi = fii(r) by demanding stability against variation of the scale: 



/J,-^-E tot (r; m, a s (n)) 
djj, 



0. 



(24) 



n=m (r) 



2. We fix the scale \x = ^{f) on the minimum of the absolute value of the last known term 
[0(a 3 s ) term] of E M (r): 



0. 



(25) 



[J,=H2 (r) 



It turns out that the total energy hardly changes whether we choose /i = /ii(r) or /x = ^(f)] 
also the scales /ii(r) and /^(r) are considerably larger than 1/r; the perturbative prediction for 
-^tot( r ) becomes unstable at r > 3 GeV -1 . 

In Fig. ^ we compare the total energy for oig(Mz) = 0.1201, 0.1181, 0.1161 with typical 
phenomenological potentials: 
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A Coulomb-plus-linear potential (Cornell potential) 



V(r) = ~ + 3 (26) 



r a 2 



with « = 0.52 and a = 2.34 GeV -1 . 
A power-law potential p6[ : 

V(r) = -8.064 GeV + (6.898 GeV)(r x 1 GeV) ai . (27) 

A logarithmic potential f2~? |: 



7(r) = -0.6635 GeV + (0.733 GeV) log(r x 1 GeV). (28) 

In order to make a clear comparison, arbitrary constants have been added to all the potentials 
and E tot (r) such that their values coincide at r = 1 GeV -1 . The total energy appears to be in 
good agreement with the phenomenological potentials in the range 0.5 GeV -1 < r < 3 GeV -1 . 
The level of agreement is consistent with the uncertainties expected from the next-to-leading 
renormalon contributions ±|Aqcd • (Aqcd^) 2 taking Aqcd = 300 MeV: indicated by the error 

bars . 

Thus, we confirm that, when the renormalon cancellation is incorporated, the QCD radiative 
corrections bend the Coulomb potential upwards at long distances. One may understand it 
similarly to the previous discussion: When the infrared cutoff ~ 1/r of Eqs. (|T2"|)-(|i"4j) is 
lowered, i.e. at large r, ats(q) grows rapidly in the integrand, which raises the total energy as 
compared to the non-running (Coulomb) case. 

Another way to understand it qualitatively is to consider the interquark force defined by 

F(r) = ~E„(r) = —V^r) (29) 

s _ CF ^m. (30) 

The last line defines the "F-scheme" coupling constant ex pip)- The interquark force is also free 
from the leading renormalon. The running of Of(/-0 is dictated by the renormalization-group 
equation: 

A^T"?"^) = Pf(<* f ), (31) 

where the first two coefficients of the beta function are universal, i.e. same as those of /3^(as)- 
When we consider effects of the QCD radiative corrections on the lowest-order Coulomb po- 
tential, one may interpret that in the QCD potential, VqcdM — —CF a s(^/ r )/ r , the coupling 
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Figure 7: A comparison of the QCD potentials calculated in ^-scheme and in F-scheme as well as the Coulomb 
potential. We set m c = and consider 4 massless flavors. The Coulomb potential is given by —Cpa/r with 



increases at long distances, so the potential will be bent downwards. This is a bad interpre- 
tation, since the QCD potential is poorly convergent at relevant distances. We should rather 
consider the interquark force. A better interpretation is that in F(r) = —CFCtF(l/r)/r 2 , the 
F-scheme coupling increases at long distances, and correspondingly |F(r)| grows with respect 
to the Coulomb force. This means that the slope of the potential becomes steeper at long 
distances. Its effect resembles an addition of a linearly rising potential to the Coulomb poten- 
tial. Thus, the effects of the radiative corrections are even qualitatively reversed, whether we 
consider Vqcd(^) or F(r) as the physically relevant quantity. 

One may verify these features in Fig. [7], in which the Coulomb potential, the V^-scheme 
potentials and the F-scheme potentials are displayed. The ^-scheme potentials are calcu- 
lated by solving the renormalization-group equation for the QCD potential. The F-scheme 
potentials are calculated by first solving the renormalization-group equation (|3lD for ap nu- 
merically and then by integrating —Fir) over r numerically; arbitrary constants are added 
such that the F-scheme potentials coincide the Coulomb potential at r = 0.4 GeV -1 . As 
can be seen, the ^/-scheme potentials become singular at fairly short distances, whereas the 
F-scheme potentials have wider ranges of validity. The 2-loop and 3-loop F-scheme poten- 
tials are consistent with the phenomenological potentials within the uncertainty expected from 
the next-to-leading renormalon contributions, in the range 0.5 GeV -1 < r < 2.8 GeV -1 and 
0.5 GeV -1 ^ r < 1.7 GeV -1 , respectively. On the other hand, the 1-loop F-scheme potential 



a = 0.279. 
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does not satisfy this criterion. 



7 Conclusions and Discussion 

For all the bottomonium states, where the predictions of perturbative QCD can be made 
reliably (i.e. as < 1), our results are consistent with the experimental data within the estimated 
uncertainties of the theoretical predictions. The obtained value for m& is in good agreement with 
the recent sum-rule calculations. The theoretical uncertainties given in Tab. || are numerically 
of the same size as Aqcd x (&xAqcd) 2 , i.e. of the effect of the next-to-leading renormalons: if 
we approximate l/a x — fJ-x, take the values of Tab. [j], and A QCD = 300 — 500 MeV, we obtain 
for the 15* state a contribution of order ±(5 — 20) MeV, for the n = 2 states a contribution of 
order ±(20 — 110) MeV and for the 3S state a contribution of order ±(50 — 250) MeV. Since the 
mass nib has been fixed on the vector ground state and has not been adjusted for higher states, 
the data at our disposal suggest that: 1) the bulk of the bottomonium spectrum is accessible 
by perturbative QCD up to some of the n = 3 states; 2) non-perturbative contributions do 
not need to be larger than 250 MeV for the reliable n = 3 states, than 100 MeV for the 
n = 2 states and than 10 MeV for m&, and they are consistent with the type associated with 
the next-to-leading renormalon. These upper bounds to the non-perturbative corrections are 
conservative and their true sizes may be considerably smaller; note that for reliable predictions 
all of \Ex P — Ex\ in Tab. |T| are smaller than 60 MeV. 

When we incorporate the cancellation of the leading renormalon contributions, the perturba- 
tive expansion of the total energy E tot {r) of a bb system, up to 0(a s ) and supplemented by the 
scale- fixing prescription (p4|) or (|25|), converges well at r <j 3 GeV -1 . Moreover, it agrees with 
the phenomenologically determined potentials in the range 0.5 GeV" 1 < r < 3 GeV" 1 within 
the uncertainty expected from the next-to-leading renormalon contributions. This establishes 
a connection between our present approach based on perturbative QCD and the conventional 
phenomenological potential-model approaches. Eventually we may merge them and further 
develop understandings of the charmonium and bottomonium systems. For instance, in the 
perturbative prediction for the bottomonium spectrum, the level splittings between the S'-wave 
and P-wave states as well as the fine splittings among the nPj states are smaller than the cor- 
responding experimental values. Although the discrepancy is still smaller than the estimated 
theoretical uncertainties of the predictions, it should certainly be clarified whether they are ex- 
plained by higher-order perturbative corrections, or, we need specific non-perturbative effects 
for describing them. On the other hand, the conventional potential-model approaches have 
been successful also in explaining the S-P splittings and the fine splittings. Hence, we expect 
that the connection would help to clarify origins of the differences of the present perturbative 
predictions and the experimental data. 

For what concerns the cc system, we obtain m c = 1243 ± 15 ± 20 ± 50 MeV from the 
mass of the J/ip state: the first error is due to the uncertainty in as(Mz), the second error is 
due to higher-order corrections, and the third error is due to non-perturbative contributions. 
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We note that this estimate is in good agreement with recent sum-rule calculations. With the 
present method, however, we cannot make reliable predictions for states higher than the ground 
state of charmonium, and, therefore, we cannot extrapolate from consistency arguments the 
size of the non-perturbative corrections. Nevertheless, from the prediction of the rj c mass, we 
may anticipate non-perturbative contributions of the order of 100 MeV to the n — 1 states. 
Also this figure is consistent with the next-to- leading renormalon effect (20 — 110 MeV). Our 
prediction for the mass of the ^(1%) state is M Bc{ls) = 6324 ± 5 ± 20 ± 40 MeV. 

We provided a novel picture on the composition of the masses of the bottomonium states. 
The picture may be contrasted to that of the QED boundstates such as positronium. For this 
system, the pole masses of a free electron and positron are well-defined. Then the boundstate 
mass is given by the sum of the pole masses minus the binding energy. For a bottomonium 
state, we cannot define the pole masses by separating b and b an infinite distance.^] Therefore, 
the only sensible description is through contributions of the gluons which reside inside the 
boundstate. In this picture: the mass of the boundstate becomes heavier than the sum of the 
quark MS masses (~ current quark masses); it is composed of the sum of the MS masses and 
the self-energies (~ constituent quark masses), supplemented by the small negative potential 
energy which binds the system together. It would be an interesting question if this picture is 
also applicable to lighter QCD boundstates. 
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